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Abstract 

Background: Chikungunya (CHIK) is currently endemic in South and Central India and exist as co-infections with 
dengue in Northern India. In 2010, New Delhi witnessed an outbreak of CHIK in the months October-December. 
This was the first incidence of a dominant CHIK outbreak in Delhi and prompted us to characterize the Delhi virus 
strains. We have also investigated the evolution of CHIK spread in India. 

Findings: Clinical samples were subjected to RT-PCR to detect CHIK viral RNA. The PCR amplified products were 
sequenced and the resulting sequences were genetically analyzed. Phylogenetic analysis based on partial 
sequences of the structural proteins E1 and E2 revealed that the viruses in the latest outbreak exhibited ECSA 
lineage. Two novel mutations, E1 K21 1 E and E2 V264A were observed in all Delhi isolates. In addition, CHIKV 
sequences from eight states in India were analyzed along with Delhi sequences to map the genetic diversity of 
CHIKV within the country. Estimates of average evolutionary divergence within states showed varying divergence 
among the sequences both within the states and between the states. We identified distinct molecular signatures of 
the different genotypes of CHIKV revealing emergence of a new signature in the New Delhi clade. Statistical 
analyses and construction of evolutionary path of the virus within the country revealed gradual spread of one 
specific strain all over the country. 

Conclusion: This study has identified unique mutations in the E1 and E2 genes and has revealed the presence of 
ancestral CHIKV population with maximum diversity circulating in Maharashtra. The study has further revealed the 
trend of CHIK spread in India since its first report in 1963 and its subsequent reappearance in 2005. 
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Findings 

Chikungunya (CHIK), a neglected vector borne disease 
until as few years ago, has grown to pandemic levels, 
with hundreds of thousands cases reported from all 
around the globe. In India, the infection reemerged in 
seven states in 2005 and by the latest report in 2010, 
has spread to more than 18 states/Union Territories 
within the country affecting more than 3.7 million 
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individuals [1]. Moreover, the number of cases has 
been grossly under reported due to mistaken diagnosis 
of dengue and non-reporting of suspected cases and 
related deaths thereby making the number of CHIK 
infection much more than reported [2]. Since its ree- 
mergence, the intensity of the infection has increased 
with every passing year with 45%-63% attack rates in 
several areas during outbreaks [3,4]. India is endemic 
to dengue fever and due to overlapping symptoms; 
CHIK infection is generally mistaken to be the former 
thereby leading to misdiagnosis. Persistent chronic 
phase is sometimes used as a distinguishing feature be- 
tween dengue and CHIK [5]. 
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Delhi witnessed a CHIK outbreak between October- 
December 2010 just after experiencing an outbreak of 
dengue. Even though CHIK has been occurring in Delhi 
since 2006 as co-infections with dengue [6] and has not 
been known to cause major outbreaks, last few months 
of year 2010 saw an unprecedented increase in number 
of affected CHIK cases, prompting us to characterize the 
strains and study the infection with special attention to 
its divergence in the country since the first outbreak in 
1963. 

Blood samples were drawn from 289 patients, who 
were clinically suspected to have CHIK, over a period of 
three months since the onset of the outbreak. The blood 
samples were received from outpatient and inpatients 
departments as part of the routine laboratory investiga- 
tions for suspected cases that were suffering from the 
classical symptoms of fever, polyarthralgia, myalgia and 
skin rashes. These patients belonged to different local- 
ities of Delhi. The sera was separated and stored at -80°C. 
All sera samples were subjected to ELISA for testing the 
presence of Chikungunya virus (CHIKV) IgM antibodies 
using the IgM capture ELISA kits supplied by National In- 
stitute of Virology (NIV), Pune to Government Hospitals 
during the outbreak. A total of 80 sera samples were used 
for the study. The sera were not subjected to passaging 
using any cell lines to ensure the sequence integrity of the 
strains. 

Viral RNA was extracted directly from sera using the 
High Pure Viral Nucleic Acid kit (Roche, Germany) 
according to manufacturer s instructions. CHIKV struc- 
tural genes El and E2 were analyzed for the present 
study. Complete gene of El (1320 bps) and E2 (1266 bps) 
and partial region (555 bps) of the El structural gene 
were amplified using new primer sets as well as previ- 
ously published primers [7,8] with modifications. For 
amplifying partial region of El gene, a nested RT-PCR 
protocol was developed to increase the sensitivity of 
RT-PCR in the clinical samples. All PCR primer sets 
used in the study and their genome positions are listed 
in Table 1. The amplicons were purified using Qiaquick 
PCR purification kit (Qiagen, USA) and sequenced using 



both forward and reverse primers. All sequences were 
trimmed, aligned and analyzed. 

Phylogenetic analyses, molecular evolution analyses of 
El and E2 genes were performed using MEGA 5 [9] and 
DNAsp softwares [10]. As references, strains of the three 
genotypes, namely ESCA, Asian and African [11] were 
taken based on their sampling dates and low passage 
history (Additional file 1: Table SI). Pair- wise genetic 
divergences within group and net evolutionary diver- 
gence between groups among the sequences were calcu- 
lated following the maximum composite likelihood 
model. The sequences were further analyzed to obtain a 
Neighbor Joining tree with bootstrap re-sampling (1000 
replicates) and validated using other distance based and 
character based tree models. Neutrality of mutation in 
the sequences was determined using Tajimas D test [12] 
and validated using Fu & Li D* & F [13]. Network ana- 
lysis was performed based on SplitsTree 4 (Ver 4.11.3) 
[14] to generate a hypothetical evolutionary pathway of 
108 Indian CHIKV samples in India since the 1963 epi- 
demic till the latest outbreak in 2010. 

Of the total 80 samples used in the study, 56 samples 
were ELISA positive and 22 samples were ELISA nega- 
tive. ELISA status of remaining two samples was not 
known. A total of 36 samples (44%) showed amplifica- 
tion of the partial gene of El. No-RT control RNA and 
RNA from healthy volunteer served as negative control. 
Of these PCR positive samples, 53% (n = 19) of samples 
were negative by ELISA, 42% (n = 15) of samples were 
positive by both methods and ELISA status of two sam- 
ples (5%) was not known. Three samples were negative 
both for ELISA and RT-PCR and were excluded from 
further analysis. Using the nested RT-PCR protocol, viral 
RNA was detected upto Day 11 from the onset of clin- 
ical symptoms (data not shown) making this an attract- 
ive alternative to the existing RT-PCR detecting 
methods. Partial region of El gene (555 bps) from 36 
samples were sequenced and aligned to study sequence 
variation. For the purpose of generating the molecular 
signature for the El gene, the complete gene was ampli- 
fied and sequenced from seven clinical isolates. A subset 



Table 1 Primer sets used in study 



Primer Name 


Sequence 5-3' 


Genome position 


Amplicon size bps 


Reference 


Where used in study 


E1-primary-F 


ACAAAACCGTCATCCCGTCTC 


10145-11158 


1013 


[7] 


Primary PCR 


E1-primary-R 


TGACTATGTGGTCOTCGGAGG 










E1-F1 


GCTCCGCGTCCmACC 


10389-10943 


555 


[8] 


Nested PCR, Phylogenetic analysis 


E1-R1 


ATGGCGACGCCCCCAAAGTC 










E1-gene-F 


GCGTACGAACACGTAACA 


9991-11310 


1320 


This study 


Molecular signatures 


E1-gene-R 


GTGCCTGCTGAACGACAC 










E2-gene-F 


AGCACCAAGGACAAOTCAAT 


8542-9807 


1266 


This study 


Phylogenetic analysis Molecular signatures 


E2-gene-R 


mAGCTGTOTGATGCAGC 
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of positive samples (n = 14) were used to amplify E2 
gene. The amplified products were purified, sequenced 
and phylogenetic analysis performed. 

Nucleotide sequence analysis was performed on partial 
regions of El and E2 genes on all Delhi strain sequences 
and compared with the 1953 Tanzanian strain (Acces- 
sion no: 045811.1). With respect to El gene, at the nu- 
cleotide level, all Delhi samples showed variations at 15 
nucleotide positions producing both synonymous and 
non-synonymous mutations (Table 2). Of the synonym- 
ous mutations, one at position C10602T was found to 
be unique to all Delhi strains. Furthermore, six strains 
showed additional variations at the nucleotide level. At 
the amino acid level, three amino acid substitutions, 
namely, El K211E, El M269V and El D284E were seen 
in all samples (Table 2). Apart from the above men- 
tioned variations, three samples showed unique amino 
acid changes, namely, IND-10-DEL48 at amino acid 



positions El V179A, El S234P, IND-10-DEL88 at pos- 
ition El R196K, IND-10-DEL108 at position El R247C 
(Table 2). Analysis of E2 gene at the nucleotide level 
revealed thirty two variations in all Delhi samples when 
compared with 1953 Tanzanian strain (Table 3). Of these 
nucleotide variations, eight resulted in amino acid changes. 
Furthermore, samples IND-10-DEL10 and IND-10-DEL15 
revealed unique amino acid variations at position E2 V50A 
and E2 C389R respectively (Table 3). 

Of the variations in both El and E2 genes, two amino 
acid residues El K211E and E2 V264A deserves special 
mention. The El K211 mutation has been reported earl- 
ier in the Asian lineage and the ECS A lineage shows El 
K211N [15]. The El K211E variation in ECSA lineage 
has been reported in one sample from Puducherry [16]. 
In the present study, all the samples collected in Delhi 
have shown this variation. Tsetsarkin and co-workers 
[17] recently studied the involvement of El -2 11 in the 



Table 2 Nucleotide and amino acid variations of El gene in Delhi samples in comparison with ECSA lineage and other 
Indian samples are shown here 


NurlpntiHp 


Amino 








Jal 1 ipiCj 


Position 


acid 

variation 


ECSA# 


Indian 
Samples 


All Delhi 
Sample 


IND-10- IND-10- IND-10- IND-10- IND-10- IND-10- IND-10- 
DEL4 DELI 7 DEL48 DEL88 DEL91 DELI 00 DELI 08 


10521 




A 






G 


10524 




C 


T/C 


T 




10529* 


V179A 


T 






C 


10545 




T 


C 


C 




10548 




C 


T/C 


T 




10580* 


R196K 


G 






A 


10602 




C 




T 


C 


10614 




G 


A/G 


A 




10623 




C 


T/C 


T 




10624* 


K211E 


A 


A/G 


G 




10641 




C 


T/C 


T 




10665 




G 






A 


10689 




A 


A/G 




G . . T 


10693* 


S234P 


T 






C 


10723 




T 


C 


C 




10732* 


R247C 


C 






T 


10734 




A 


C/A 


C 




10743 




G 


A/G 


A 




10746 




A 


G 


G 




10798* 


M269V 


A 


G/A 


G 




10845* 


D284E 


T 


A/T/C 


A 




10848 




G 






A 


10893 




A 


G 


G 





Asterisk next to nucleotide position denotes non-synonymous mutations. Dot represents same nucleotide as ECSA. Blank represents same nucleotide as All Delhi 
samples. 

#ECSA strains (Accession nos. AF369024.2 and HM04581 1.1). 
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Table 3 Nucleotide and amino acid variations of E2 gene 
in Delhi samples in comparison with ECSA lineage and 
other Indian samples are shown here 

Samples 

ECSA# Indian All Delhi IND-10- IND-10- 
Samples Samples DEL10 DELI 5 



Nucleotide Amino 
Position acid 

variation 



8690* 


V50A 


j 






8703 




j 


Q 


Q 


8721 




j 


C/T 


Q 


8727 




j 


C/T 


Q 


8748 




Q 


T/C 


J 


8757 




J 


C/T 


Q 


8784 




Q 


G/T 


Q 


881 1 






Q 


Q 


8814 




c 


T/C 


J 


8859 




Q 


A/G 




8865 






G/A 


Q 


8910 




c 


T/C 


J 


8952 






T/C 


J 


9031* 


A164T 


Q 


A/G 




9081 




Q 


A/G 




9138 




J 






9173* 


121 IT 


J 






9270 




Q 






9276 




Q 


j 


J 


9285 




Q 


T/C 


T 


9332* 


V264A 


J 


J 


Q 


9437* 


S299N 


Q 






9476* 


T312M 


Q 


T/C 


T 


9480 




C 


T 


T 


9489 




G 


A/G 


A 


9525 




T 


C 


C 


9571* 


A344T 


G 


A 


A 


9588 




C 


T/C 


T 


9609 




C 


T 


T 


9610 




T 


C/T 


C 


9615 




C 


T 


T 


9664* 


S375T 


T 


A/T 


A 


9698* 


V386A 


T 


C/T 


C 


9706 


C389R 


T 






9741 




C 


T/C 


T 



Asterisk next to nucleotide position denotes non-synonymous mutations. Dot 
represents same nucleotide as ECSA. Blank represents same nucleotide as All 
Delhi samples. 

adaptation of the Asian strain and eventually in its 
spread. They illustrated that El -2 11 along with another 
mutation at position El -98 were associated with CHIKV 
sensitivity to Ae. albopictus-zd&ptive effect of El A226V 



mutation. Similarly, in E2, amino acid variation E2 
V264A was seen in all samples of Delhi population. Both 
the El K211E and E2 V264A variations have also been 
reported in the autochthonous imported case in France 
from India during the same period in 2010 from a differ- 
ent region [18] as well as from cases in Tamil Nadu and 
Andhra Pradesh [19]. 

Furthermore, critical evaluation of the complete El 
and E2 gene sequences data allowed us to generate a 
molecular signature for the specific genotypes of CHIKV. 
The molecular signature illustrate that the Delhi samples 
showed a different signature from those seen in the 
reference strains and other strains from India reported 
since 1963 and the latest documented epidemic in 2005 
onwards (Tables 4 & 5). A recent study has reported 
emergence of a novel subgroup among CHIKV strains 
in Asian countries; however, homologous recombin- 
ation did not contribute to the genetic diversity of 
these strains [20]. It should be noted that the above 
study did not include any of the strains from the latest 
outbreak in 2010 in India. Molecular signatures identi- 
fied by our study raises a possibility of recent recom- 
bination events in the Indian CHIKV population and 
warrant further detailed analysis. 

Phylogenetic analysis based on El of all Delhi sam- 
ples with strains previously reported from eight states 
in India since 1963 (n = 109) revealed that all Delhi 
samples were of ECSA lineage and formed a single 
subclade within the Indian sequences (Figure la). Simi- 
lar result was seen with E2 gene (n = 42) (Figure lb). 
Evolutionary Divergence over 102 sequence pairs for 
El gene was calculated with a total of 143 positions in 
the final dataset and sequences from Maharashtra were 
found to be most divergent among all sequences ana- 
lyzed. Tests of neutrality determined by Tajimas D test 
showed positive selection in El gene, that was further 
validated by Fu & Li's D* and F* tests (p < 0.05). Previ- 
ous studies [15] have shown that El K211 is a posi- 
tively selected site with a posterior probability of 
>75%. Another very recent study conducted on sam- 
ples collected from different regions in south India, 
also showed that El K211E mutation is being positively 
selected [19]. With respect to fitness of these strains in 
the vector, it should be noted that Delhi had tradition- 
ally predominance of Ae. aegypti population; however, 
it has been shown recently that Ae. albopictus is 
adapted to urban breeding conditions and also plays a 
role in dengue transmission in Delhi and National Cap- 
ital Territory of Delhi (NCT) [21], It is noteworthy that 
Ae. aegypti is very prevalent in Tamil Nadau and 
Andhra Pradesh from where the samples in Sumathy 
et aL, study [19] were taken while Ae. albopictus dens- 
ity is high in the neighbouring states. Since it has been 
shown in previous studies [22] that mutations in E2 
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Table 4 Molecular signatures derived on amino acid residue variations in chikungunya whole El proteins are shown 
below 



Collection Year Sample El protein: Amino Acid Residue Position 







72 


98 


142 


145 


162 


211 


225 


226 


269 


276 


284 


296 


315 


322 


1966/2005 


Africa 


S 


A 


1 


A 


V 


K 


A 


A 


l/V 


1 


D 


V 


V 


A 


1953 


ECSA 


N 


A 


1 


T 


1 


K 


A 


A 


M 


M 


D 


L 


V 


A/V 


1985/2007 


Asia 


S 


T 


l/V 


A/S 


1 


E 


S 


A 


M 


M 


D 


L 


V 


A 


2010 


Delhi* 


N 


A 




T 




E 


A 


A 


V 


M 


E 




V 


A 


1963 


West Bengal 


S 


T 


V 


S 




E 


S 


A 


M 


M 


D 




A 


A 


2003/2006 


Maharashtra 


N/S 


A/T 


l/V 


T/S 




K/E 


A/S 


A 


V/M 


M 


E/D 




V/A 


A 


2006 


Gujarat 


N 


A 




T 




K/N 


A 


A 


V 


M 


E 




V 


A 


2006 


Puducherry 


N 


A 




T 




K/E 


A 


A 


V 


M 


E 




V 


A 


2006/2007/2008 


Kerala 


N 


A 




T 




K 


A 


A/V 


V 


M 


E 




V 


A 


2007 


Uttar Pradesh 


N 


A 




T 




K 


A 


A 


V 


M 


E 




V 


A 


2006/2008 


Karnataka 




A 




T 




K 


A 


A/V 


V 


M 


E 




V 


A 


2008/2009 


Andhra Pradesh 


N 


A 




T 




K/N 


A 


A 


V 


M 


E 




V 


A 



Amino acid residue positions of the complete E1 proteins, the Country/State and year of collection of CHIKV are denoted. (-) denotes sequence information not 
available. 



play an epistatic role on El gene, it will be very inter- 
esting to analyze the effect of these mutations with re- 
spect to vector adaptability and infectivity in the presence 
and absence of El A226V mutation which also is posi- 
tively selected [19]. Estimates of average evolutionary 
divergence (AED) over sequence pairs within the states 
showed varying divergence among the sequences both 
within the states and between the states. Sequences 
derived from samples from Karnataka and Uttar Pradesh 
showed least divergence, while Maharashtra showed max- 
imum divergence suggesting existence of ancestral popu- 
lation in Maharashtra. 



Following phylogenetic analysis and genetic diversity 
analyses, attempts were taken to construct a network on 
the basis of the important mutational events of the virus 
along its geographical spread along with the sampling 
date. The purpose of the network was to visualize the mu- 
tational paths based on the variations in the sequence ana- 
lyzed as has been ascertained in recent studies [23,24]. In 
the present study, the mutation status combined with 
sampling date and geographical information provides 
insight in the spread of the strains with the country since 
the first incidence of the infection. Since the software also 
predicts hypothetical ancestral variants, this tool is useful 



Table 5 Molecular signatures derived on amino acid residue variations in chikungunya whole E2 proteins are shown 
below 



Collection Year Sample 



E2 protein: Amino Acid Residue Position 







2 


118 


149 


157 


205 


211 


255 


264 


299 


317 


318 


344 


375 


384 


386 


1966/2005 


Africa 


T 


S 


K 


V 


G 


T 




V 


N 


V 


T 


T 


S 


T 


V 


1953 


ECSA 


T 


S 


K 


V 


G 






V 


S 


V 


V 


A 


S 


M 


V 


1985/2007 


Asia 




G 


R 


A 


D 


T 


V 


V 


N 




R 


T 


s 


V 


V 


2010 


Delhi* 


T 


S 


K 


V 


G 


T 




A 


N 


V 


V 


T 


T 


M 


A 


1963 


West Bengal 




G 


R 


A 


D 


T 




V 


N 


V 


R 


T 


s 


V 


V 


2003/2006 


Maharashtra 


l/T 


G/S 


K 


V 


D/G 


T 




V 


N 


V 


R/V 


T 


S/T 


V/M 


V/A 


2006 


Gujarat 


T 


S 


K 


V 


G 


T 




V 


N 


V 


V 


T 


T 


M 


A 


2006/2007/2009 


Kerala 


T 


s 


K 


V 


G 


T 




V 


N 


V 


V 


T 


T 


M 


A 


2006 


Rajasthan 


T 


s 


K 


V 


G 


T 




V 


N 


V 


V 


T 


T 


M 


A 


2006 


Tamil Nadu 


T 


s 


K 


V 


G 


T 




V 


N 


V 


V 


T 


T 


M 


A 


2006 


Karnataka 


T 


s 


K 


V 


G 


T 




V 


N 


V 


V 


T 


T 


M 


A 


2006/2008/2009 


Andhra Pradesh 


T 


s 


K 


V 


G 


T 




V 


N 


V 


V 


T 


T 


M 


A 



Amino acid residue positions of the complete E2 proteins, the Country/State and year of collection of CHIKV are denoted. (-) denotes sequence information not 
available. 
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GU562826.1|Kerala|2009 

GU562829.1|Kerala|2009 

GU562828.1|Kerala|2009 

GU562830.1|Kerala|2009 

GU562827.1|Kerala|2009 

EU372006.1|Kerala|2007 

GQ428210.1|Kerala|2006 
EF027134.1|Andhra Pradesh|2006 
EF027137.1 |Rajasthan|2006 
HM045794.1 |Gujarat|2006 
JF272479.1|Gujarat|2006 
EF027136.1 |Maharashtra|2006 
EF027138.1|Tamil Nadu|2006 
JF274082.1|Gujarat|2006 
FJ000066.1|Kerala|2006 

♦ lND-10-DEL75|2010 

♦ lND-10-DEL7|2010 

♦ lND-10-DEL5|2010 

♦ lND-10-DEL84|2010 

♦ lND-10-DEL91|2010 

♦ lND-10-DEL4|2010 

♦ lND-10-DEL12|2010 

♦ lND-10-DEL100|2010 

♦ lND-10-DEL9|2010 

♦ lND-10-DEL102|2010 

♦ IND-10-DEL1 12010 

♦ lND-10-DEL123|2010 

♦ lND-10-DEL10|2010 

♦ lND-10-DEL15|2010 
HQ529779.1|Andhra Pradesh|2009 
HQ529777.1|Andhra Pradesh|2008 

GQ428211.1|Kerala|2006 
EF027135.1|Karnataka|2006 
HM045811.1|Tanzania|1953 
- AF369024.2|S27 
HM045813.1|West Bengal|1963 
AY424803.1|Maharashtra|2003 

EU192143.1|lndonesia|2007 
HM045797.1|RSU|1985 



HM045817.1|Senegal|2005 
HM045816.1|Senegal|1966 



M20303.1|ONyong-nyong virus 



Figure 1 A & B: Phylogenetic relationship among CHIKV isolates from eight states of India. A) based on partial E1 gene sequences (428 nt) 
performed on 109 nucleotide sequences and B) on partial E2 gene sequences (1002 nt) performed on 42 nucleotide sequences. The tree was 
constructed using Neighbor-Joining method. The percentage bootstrap support values (1000 replicates) are shown next to the branches. The tree is 
drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary 
distances were computed using the Kimura 2-parameter method and are in the units of the number of base substitutions per site. Evolutionary 
analyses were conducted in MEGA5. All sequences of strains sequenced for this study has been submitted in GenBank; other sequences were retrieved 
from GenBank. Scale bar indicates number of base substitutions per site. Onyong-nyong virus used as out-group. Isolates sequenced in this study are 
indicated by 



to predict possible mutational steps in the absence of sam- 
pling. The network created using El gene sequences 
(n = 108) from different states of India shows evidence of 
possibly an independent introduction of CHIKV strains of 
Asian and ECSA lineage into the country with the Asian 
lineage introduced first followed by the ECSA lineage 
(Figure 2) also been discussed in previous reports [25]. 
The network also clearly shows Delhi strains forming a 
separate sub-clade along with one strain from Pudu- 
cherry thereby indicating population amplification of 
the virus from Puducherry in Delhi. All our analyses 



were performed on partial El gene due to the availability 
of maximum number of sequences in this gene region 
(n = 108). To validate our result, we generated the network 
on complete gene of El gene with six strains from Delhi 
and 38 strains from other parts of India. The results cor- 
roborated with the network derived from the partial El 
gene sequences (data not shown). 

In summary, our study based on the recent outbreak 
in Delhi has revealed important insights to the nature of 
the spread of CHIKV within the country. Even though 
CHIKV of the latest outbreak displayed ECSA lineage, 
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Figure 2 Pruned quasi-median network of CHIKV partial El gene sequences. 428 nucleotides of E1 gene (n = 108) was used to generate the 
network, here a specific color code are assigned to each state/Nation and diameter of nodes corresponds to the sample size within each node. 
The edges joining the nodes are not in scale. Hypothetical ancestral strains or strains not sampled are also represented as small nodes white in 
color. The proposed evolutionary path is shown by red arrows. Sky blue star represents substitution E1 K21 1E. 
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there is emergence of a distinct molecular signature 
within the strains emphasizing that accumulation of 
mutations in the viral genome is leading to appearance 
of new subgroups and suggests a dynamic evolution of 
the virus [18]. Understanding the consequences of such 
changes in viral genomes is vital to prepare for any pub- 
lic health disaster like a serotype change or more serious 
clinical implications. Our attempts to trace the path of 
Indian CHIKV has also revealed that ECSA lineage may 
have emerged in India much later than the Asian lineage 
and have flared up along with the Reunion Island pan- 
demic. Our study further shows evidence that important 
mutational events are also taking place within this clade 
that may have a potential role in the evolution of the 
virus and necessitates further detailed studies. 
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